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I Pion elastic scattering on deuterium is studied in the KMT multiple scatter- 

ing approach developed in momentum space. Using a Paris wave function and 
the same methods and approximations as commonly used in pion scattering 



on heavier nuclei excellent agreement with differential cross section data is 



obtained for a wide range of pion energies. Only for T w > 250 MeV and very 

H : 

backward angles, discrepancies appear that are reminiscent of disagreements 
in pion scattering on 3 He, 3 H, and 4 He. At low energies the second order cor- 
rections have been included. Polarization observables are studied in detail. 
While tensor analyzing powers are well reproduced, vector analyzing powers 

exhibit dramatic discrepancies. 
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I. INTRODUCTION 



In the realm of pion-nuclear physics, the pion interaction with the deuteron is of special 
interest 0] . On the theoretical side, it is the simplest few-body system in which conventional 
theories (Faddeev theory, coupled ttNN — NN equations and multiple scattering theory) can 
be rigorously tested. The knowledge obtained from these studies can then be extended to 
many-body systems. On the experimental side, a large set of measurements for differential 
cross sections and polarization observables are available in red scattering. Thus, there is 
sufficient motivation to study pion interaction mechanisms with very light nuclei in detail. 

In recent years the three-body approach has been used extensively in describing the 
pion-deuteron interaction. From the point of view of solving the ttNN three-body problem 
directly this is clearly an important achievement. However, due to the complexity of the 
three-body formalism the application of the obtained results to the case of pion interaction 
with heavier nuclei may be more difficult. 

On the other hand, a microscopic description of the pion-nuclear interaction in the frame- 
work of the multiple scattering theory p],[3| has continuously been improved. More recent 
developments are due especially to techniques developed to perform calculations in mo- 
mentum space using the KMT || formulation of the multiple scattering theory. In mo- 
mentum space nonlocalities of the pion-nuclear interaction, off-shell extrapolations of the 
pion-nucleon scattering amplitudes and exact treatment of Fermi motion can be taken into 
account The KMT approach is important to avoid double counting of pion rescatter- 

ing on the same nucleon. With the addition of the phenomenological p 2 term |§, which is 
responsible for real pion absorption and second order effects, the momentum space formal- 
ism was successful in describing the pion-nuclear interaction in the A resonance region as 
well as at low energies for a large set of nuclei with A=4-40. In Refs. 0, this method has 
been extended successfully to the description of pionic atoms as well. 

The aim of the present work is a systematic investigation of pion scattering on the 
deuteron in the region of pion kinetic energies of T n =60-300 MeV. We employ a multiple 
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scattering framework in momentum space using only elementary amplitudes extracted from 
ttN scattering data and realistic deuteron wave functions. Due to the latter ingredient we 
are in a good position to fix the nuclear structure input in order to shed more light on the 
reaction mechanism. Another purpose is to develop a pion deuteron interaction that can be 
used for pion photoproduction on the deuteron in a straightforward way. Such an approach, 
applied to pion scattering and pion photoproduction on the trinucleon [|7||| has provided an 
excellent description of the differential cross section. This has motivated us to extend our 
momentum space approach for these reactions to the deuteron. 

Our study has been divided into two parts. In this paper we investigate pion elastic 
scattering on the deuteron, testing our multiple scattering approach by comparing with all 
available experimental data. Coherent tt° photoproduction is considered in the second part. 
The main aspects of our formalism based on the KMT multiple scattering approach and 
coupled-channels method are given in section 2. Section 3 presents our results for elastic 
pion-deuteron scattering while our conclusions are summarized in section 4. 

II. FORMALISM 

A. General overview 

In multiple scattering theory the pion-nuclear T-matrix can be presented as 

T = J2T l (E), (1) 

i=l 

where 

T^E) = n (E) + Ti(E)G (E) ^Tj(E) (2) 

and Ti(E) is defined solution of the equation 

t(E)=v + vG (E)t(E). (3) 

In Eqs. (2-3), Gq{E) is the Green's function of the non- interacting pion-nuclear system 
and v is the potential of the pion-nucleon interaction. Note that the r-matrix in Eq. (3) 
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differs from the standard two-body t-matrix for free pion-nucleon scattering through the 
Green's function Gq{E) which contains the many-body nuclear Hamiltonian. For the free 
pion-nucleon t-matrix we have the equation 

t{u) = v + vg (uj)t(uj) (4) 

with the free pion-nucleon Green's function go(u>). 

Following the KMT- version of multiple scattering theory ||, Eqs. (1-4) are equivalent 
to the system of integral equations for the auxiliary matrix T' = ( A ~^ T 

T\E) = U\E) + U'(E)G (E)PT'(E) , (5a) 
U'(E) = (A - l)r(E) + (A - l)r(E)G (E)QU(E) , (5b) 

r(E) = t(u) + t(u)[G (E) - 5 oH](P + Q)t(E) , (5c) 

where the operators P =\ >< | and Q = J2 m ^o \ m >< m \ project the nuclear state 
vectors into the ground state | > and all possible excited states | m >, respectively In the 
case of the deuteron, the operator Q projects onto the states of the continuum. 

One of the main problems which arise in applying these equations to pion deuteron 
scattering is the correct description of the coupling to the break-up or continuum channels. 
In this case the Lippmann-Schwinger equation has a noncompact kernel and, therefore, 
cannot be solved numerically. It is well known that the solution to this problem was found 
by Faddeev who from the system of Eqs. (5) derived a new set of equations which do not 
have this inherent shortcoming. 

In this paper we demonstrate that multiple scattering theory is able mainly to reproduce 
existing experimental data using several standard approximations for the solution of the 
Eqs. (5). This in turn generates a simple pion-deuteron interaction that will be used in 
our description of coherent tt° photoproduction, discussed in the following part. As a first 
step, we neglect the contribution from the coupling to the break-up channel. This means that 
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in Eqs. (5b-5c) the contribution from the Q-operator is dropped. The motivation for this 
approximation is that contributions from the diagonal matrix elements dominate due to the 
large overlap of the nuclear wave functions in the initial and final states. As we will see below 
this approximation is a resonable starting point for the study of pion-nuclear interaction. 
However, at low energies where for isoscalar nuclei the contribution from the P operator is 
very small, the coupling to the break-up channels becomes important. 

The next approximation is the special choice for connecting the pion-nucleon energy, 
u, with the pion-nuclear one, E. In general, the energy of the pion-nucleon system is a 
dynamical variable in the many-body system [K|. Thus, there are many possibilities for 
relation uo and E; a detailed investigation of this issue was done in Ref. |J. Here we follow 
the prescription of Landau and Thomas |4]] which is called modified impulse approximation 
with the three-body choice for the reaction energy. The relativistic generalization of this 
approximation in the case of the deuteron leads to the following expression for the pion- 
nucleon energy || 

T' + q) 2 



U) = E + 171k + Mn — 



(6) 
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where and are the pion and nucleon masses, respectively. As shown in Refs. PJ1T 
this choice minimizes the contribution from the second term in Eq. (5c). 

The scattering amplitude F MfM .(q' ,q) is connected with the T'-matrix in Eq. (5a) by 
the relation 



M( q')M(q) A 
2tt A - I 



F MfMi (q',q) = - 1 \; — - A , < 7r(q'),f | T\E) | i,7r(g) >, (7) 



where q and q' are the pion momenta in the intial and final states, respectively, | i >=\ 
1 + Mj > and | / >=| 1 + Mf > denote the nuclear initial and final states, with the nuclear 
spin projection Mj in the initial and Mf final states. The pion-nuclear reduced mass is 
given by M.(q) = E n (q)EA_(q)/E(q), where E(q) = E n (q) + EA_(q) is the total pion-nuclear 
energy. To shorten the notation in Eq. (7) we skipped the dependence on the energy E for 
the scattering amplitude. For a given pion energy, E is fixed and in the following we use 



E = E(q). In the same way we can express the pion-nuclear potential in momentum space 
VM f Mi(q',q) via matrix U'{E) from Eq. (5b) Then, in accordance with Eq. (5a) the elastic 
scattering amplitude can be constructed by solving the integral equation with relativistic 
kinematics 

F MfMi (<l .9) " ,9) " ^ E(g)-W)+^ ' ( } 

where the factor a = (A — 1)/A is important to avoid double counting of pion rescattering 
on the same nucleon, which is already included in the free pion-nucleon t-matrix. Note that 
due to this factor, Eq. (8) is not the standard Lippmann-Schwinger equation. 

In accordance with the approximations discussed above and after neglecting the contri- 
butions from the second terms in Eqs. (5b) and (5c) the momentum space potential of the 
pion-nuclear interaction can be given by 

V M fMi(q', q) = Vcoui.(<f - T, R)S Mf M t + V$ fM .(q', q) , (9) 

which contains the Coulomb potential in momentum space, cut-off at radius R. The so-called 
first-order potential of the strong pion-nuclear interaction V^ /M . is related to the free ttN 
scattering t-matrix 



<W.g) = J M{q ^ M{q) < n(q')J | J2 *» \h^q)>, (10) 
The free pion-nucleon scattering t-matrix is defined in the following way 

Un =< q',p' | t(u) \p,q>= 1 = f wN (uj,9*)5(p' + q' -p- q) , (11a) 

^H{q',p')li(q,p) 

f nN {uj, 01) = A + A T t- r + ia-fa x $ f ](A s + A ST t- r) , (lib) 

where fj,(q,p) = E^{q)E^{p)/uj is the pion-nucleon reduced mass, p and p' are the nucleon 
momenta in the initial and final states (in pion-nuclear cm. system). In the pion-nucleon 
scattering amplitude f n ^, which is in general relativistically invariant, q { and are the unit 
vectors for the initial and final pion momenta in the ttN cm. system. The vectors f and t 
are the usual isospin operators for the target nucleon and pion. 
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The scalar functions Ai(u, cos6**) (i = 0,S,T, and ST), that depend on the total pion- 
nucleon energy cu and the pion angle 6* in the irN cm. system, are the usual combinations 
of partial rrN scattering amplitudes fj- (u) and Legendre polynomials Pi n (cos9*), where 
is the pion-nucleon angular momentum. The definition of these amplitudes has been given in 
our previous work |]]. Note that their on-shell value was obtained from phase-shift analysis 
of the elastic pion-nucleon scattering data. 

The off-shell extrapolation of the ttN partial amplitudes was constructed using a sepa- 
rable form 



ft\oS- shell) -ft^f-^Qf^ 1 (12) 



with ttN form factors 

^ - irrhr ■ (13) 

Here we employ the value r = 0.47/m, consistent with the analysis of the separable ttN 
potential of Ref. JTJ 

Finally, we mention a useful approximation related to the treatment of nucleon Fermi 
motion. In accordance with the results of Refs. [P|,|i"5| the internal nucleon momentum can 
be replaced by an effective value 

q 1 

P^Peff. = -^-^(q-q')- (14) 

This choice is based on the fact that in the case of Gaussian wave functions, which can 
reproduce the dominant S-wave part of the deuteron ground state at low momenta, such a 
replacement treats the linear p/2M terms in the pion-nucleon scattering amplitude exactly. 

B. Partial wave decomposition and polarization observables. 

We express the scattering amplitude in terms of partial amplitudes using the represen- 
tations of total angular momentum j with projection m 
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FM fMi {q',q) = ^(2j + l)i L ^Y^{q>)F^ ^{q> , q)Y^(q) 



K 1 j\(L n l j 
Ml Mf -mj \Mn Mi -m 



(15) 



where the Y^f [q) are the spherical harmonics for the pion waves and L w (Z^.) is the an- 
gular momentum of the incoming (outgoing) pions. Note that due to parity conservation 
(— l) L ' r+L ' r = 1. For the amplitude VmsM, we perform an expansion identical to Eq. (15). 

Substituting the above expansions of Vm^m* and F MfM . into Eq. (8) we obtain the 
following system of integral equations for the partial wave amplitudes 

n, L M,*) = Vl^(q,q)--Zj M(qf/) * E{q) _ E(qlf)+ie dq . (16) 

This equation is solved using the matrix inversion method. In our evaluation of the par- 
tial amplitudes we also take into account the Coulomb interaction applying the matching 



procedure developed by Vincent and Phatak ||14j| . 



Expressions for the polarization observables can be obtained using the spherical tensor 
operator r^. K which is defined as 

<jM>\^\jM> = rk(-iy^( J M _ J M , k ), (it) 



were we use the notation J = \JU + 1. Analyzing powers can then be expressed as expec- 
tation values 

Tr(Fr kK Fl) 

hK ~ Tr(FFt) ' 1 } 

There are four polarization observables which can be measured in nd elastic scattering 
using a pion beam and a polarized deuteron target. These are the vector analyzing power 
iT u , and the three tensor analyzing powers T 20 , T 21 and T 22 . The last two can be measured 
only in linear combinations, such as 

T21 =T 21 + ^(T 22 + T 20 /V6), (19a) 
t 22 = T 22 +T 20 /V6. (19b) 



The differential cross section and four polarization observables can be expressed in terms 
of five amplitudes, called A, B, C, V and £ , which correspond to our amplitudes F MfM . with 
different initial and final deuteron spin projections. Following Robson's convention |15j they 
are defined as 



'M f Mi 



F++ F 0+ F_ + 
F+o F 00 F_ 
F+- -F _ F — 



ABC 
V £ -V 
C -B A 



(20) 



where the sign +, 0, — corresponds to the deuteron's spin projection = +1, 0, —1. Then 
the five observables can be written as 



da 2g' /: 2 



dQ 3q 



A 



.6 | 2 + | C | 2 + | V | 2 +- | £ | 2 ) = ^a(^) 

2 



(21a) 



iT u = J-Im[V*(A -C)- B*£) }/a{6 n ) 



(21b) 



T 20 = ^(| A | 2 + | B | 2 + | C | 2 -2 | V | 2 - | £ | 2 )/a(^) 



(21c) 



T, 



21 



'-Re[V* (A-C) + B*£) ]/a{6 n ) 



(21d) 



T 22 = -j=[2Re{A*C)- \B\ 2 ]/a(9 n ) 



Note that due to the relation 15 



£ = (A - C) - V2(B + V) cot 9 n 



(21e) 



(22) 



only four amplitudes are independent. 

In our numerical calculations we will use a system of coordinates which corresponds to 
the Madison Convention: A right-handed coordinate system in which the positive z-axis is 
along the beam direction (along the initial pion momentum q) and the y-axis is along the 
vector [q x q']. 
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C. Pion- nuclear potential 



Here we discuss the first-order potential vj^j M ., defined in Eq. (10). It can be rewritten 



as 



V^ fMi (q',q) = -2 



^M{q')M{q) 



J dr^* Mf i 



r) cxp i 



,iQ ■ r, 

2tt J "' ~ lvl f y ~ ' ' y 2 ' 

where Q = q — q' is the momentum transfer. 

The deuteron wave function in this expression is given by 



t wN (uj;q',q)^ Mi (r), (23) 



(24) 



where I = and 2 denote the S- and /^-components, respectively, xtn s denotes the spin 
wave function for the two-nucleon system with total spin s = 1 and projection m s . After 
a multipole decomposition for the exponential, exp(^ : ), in Eq. (23) the pion-nuclear 
potential can be expressed via reduced nuclear matrix elements M SLJ {Q) defined as 



</IE {y^)®* 

3=1 



J . ,Qr, 



m Jl( - 2 



)\i >= ("I) 



l-M, 



1 J 1 

M f M Mi 



M SL j(Q) , (25) 



where ]l{.z) is the spherical Bessel function, a s = 1 and a for S = and 1, respectively. 
Using standard techniques one finds 



MslAQ) = 3\/^7 J -J2U'LJ 

v 47r VI 







1 V 1 
1 I 1 
J L S 



) R ( in\Q)- 



(26) 



The radial integrals R^(Q) are defined as 



i#>(Q)= J drU v {r) 3L {9L )Ul{r) 



(27) 



and satisfy the normalization condition 



R ( S(Q) + R { S(Q) Q =i 



(28) 



Explicit expressions for the matrix elements Mslj(Q) are given in the Appendix. 
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The radial wave functions Ui(r) are parameterized |16j as a discrete superposition of 
Yukava-type terms 



N 



u o(r) = J2 C n ex P (~ m n r) , 



n=l 



N 



3 3 

U 2 (r) = ^ D n exp(-m n r)(l + —— + 

n=l 



m„r mlr 2 ' 



( 3 Si — component) (29a) 



^Di — component) (29b) 



where the 13 coefficients C n , D n and masses m n have been calculated following the prescrip- 



tion given in Ref. [16 



In order to facilitate a qualitative understanding of the behavior of the polarization 
observables we present expressions for the five amplitudes A,...,£ in the plane wave ap- 
proximation (PWIA) (the second term in Eq.(8) is dropped). Taking into account only 
contributions from the S'-wave component and their interference with the .D-wave compo- 
nent of the deuteron wave function (neglecting terms with V — I — 2), we obtain 



A = 2 A W 



It 



(0) 
00 



B 



sin 6L 



(1—3 cos 9 n ) R ( 



(2) 
02 



1 



(2) 



V2A S W S lR$ + —R 02 



3A W R$ 



W-(l + cos^) A W R { 



(2) 
02 5 



V = smO n V2A S W S (b$ + ±=B$ J + 3 A W R t 



(2) 
02 



£ = 2 A W 



J$ + -L (1 - 30080*) Jig 



(30a) 
(30b) 

(30c) 
(30d) 
(30e) 



Here, and are the isoscalar scalar and spin-flip 7riV-scattering amplitudes from Eq. 
(11), and Ryi are the radial integrals defined in Eq. (27), the coefficients Wo and Ws are 
kinematical factors arising from the Lorentz transformation of the elementary amplitude 
from the irN cm. to the nd cm. system. They can easily be obtained from Eqs. (10- 
11) and the angle transformation of [<& x qf] in the elementary amplitude Q]. The complete 
expressions of the five amplitudes (including the full contribution from the deuteron D-state) 
are given in the Appendix. 
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III. RESULTS AND DISCUSSION 



A. Pion scattering on unpolarized targets 

We begin our discussion with some of the main features of the pion-nuclear interaction. 
The results of our calculations for pion kinetic energies in the lab system, T 7r = 65, 181 and 
254 MeV, are shown in Fig. 1. One of the most important properties of the 7rN interaction 
in this energy region is the dominance of the p-wave contribution coming from the A-isobar 
excitation. This feature is reflected in the coherent scattering process which is proportional 
to A (nuclear mass number) and described by the scalar-isoscalar part {Aq) of the itN 
amplitude. Since thep-wave part of this amplitude has a cos 6^ dependence, the differential 
cross section experiences a minimum around 8 n = 90° (see dashed curves). The position of 
this minimum is slightly shifted from the Lorentz transformation of the pion angle from the 
ttN to the ird cm. frame and from the s- and p-wave interference. The spin-flip transition 
coming from the amplitude As (which is also mainly of p-wave nature), is proportional to 
sin 9^. The corresponding spin-flip contribution fills in the minimum, as shown in Fig. 1. 

In the A-resonance region (2^-= 180 MeV) our results are in fairly good agreement with 
the experimental data and are basically the same as the results of three-body Faddeev 
calculations [|IJ . This may be due to the dominance of the A-resonance contribution. It had 
been shown in Ref. |T7[ that in this case the two-body approach is a good approximation 



for the more elaborate three-body Faddeev framework. 

However, at lower energies (T v = 65 MeV) our model fails to reproduce the experimental 
data in the forward direction. This is related to the well-known problem of describing the 
low-energy S'-wave pion-nucleon interaction in nuclei with zero isospin. In this case, the 
contribution from the large isovector ttN scattering amplitude At cancels and only the 
small Aq amplitude remains in the first-order potential. Therefore, in the zero-energy limit 
our approach gives a small value for the 7T(i-scattering length 

a$ = lim « -0.015/m , (31) 
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where 5^ is the S'-wave red scattering phase shift calculated with the first-order potential. 



The experimental value is about five times larger: a^f' = —(0.073 ± 0.02)/m |18| . In such 
a situation one expects higher-order effects in the pion-nuclear interaction to be important. 



The well-known result (see for example Refs. [19,20]) for the second-order correction is 



a% = 2(Al-2A* T )Q) C \ (32) 

where < 1/r >= 0.46/m" 1 and C = (1 + m 7r /M)(l + m 7T /2M)~ 1 . In this expression the 
contribution from the isovector amplitude At which describes double isovector scattering 
(including charge exchange) becomes important. Taking into account such a correction for 
the scattering length we obtain 

a nd « -(0.015 + 0.034)/m = -0.049/m (33) 

which is in better agreement with the experimental value as well as with results from Faddeev 
calculations: a^ d ' = 0.046/m |HJ (for the same set of A and At amplitudes). As shown 
in Fig. 1, this correction factor also improves the agreement with the experimental data for 
the differential cross section at 7^=65 MeV (see dash-dotted curve). 

The role of the 3 /^-configuration in the deuteron wave function is also illustrated in Fig. 
1. Its contribution can be seen by turning off the nuclear matrix elements of Eq. (26) with 
L = 2 (see dotted curves in Fig. 1), thus retaining the normalization condition of Eq. (28). 
As expected, the .D-state contribution becomes significant at large angles, 9 n > 90°, for pion 
energies T n >180 MeV which correspond to momentum transfers of Q > 1.78/m^ 1 . 

Fig. 2 presents our results for higher pion energies. Note that with increasing pion 
energies the second-order correction discussed above becomes less important due to the 
dominance of the pion-nucleon p-wave contribution. However, it is still visible at T 7r =142 
MeV and it improves the description in the forward direction (the corresponding result is 
not shown in Fig. 2). Fig. 2 illustrates the role of pion rescattering contributions coming 
from the second term of Eq. (8) by comparing the plane wave (without second term of 
Eq. (8)) calculations (dashed curves) with the full results (solid curves). Pion rescattering 
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becomes very important in the A-region (140 < T n < 250 MeV), especially at large angles. 
Here it can reduce the plane wave results by a factor of two. However, with further increase 
of energy the 7riV-interaction becomes weaker again, and, consequently, the rescattering 
contributions become smaller. 

At T n > 250 MeV the main problem which arises in our approach (as well as Faddeev 
calculations) is the discrepancy with experimental data at backward angles. In contrast to 
Faddeev results which usually overestimate the data in this region by a factor of two [[!]], 
our calculations tend to underestimate the measurements. We point out that a similar 
disagreement has been found for pion scattering on 3 He and 4 He |6|,^2|. 

In the framework of Faddeev equations there have been numerous attempts to improve 
the situation at backward angles, mainly by changing the two-body inputs. One of the 
more controversial attempts addressed the Pn part of the ttN amplitude with its division 
into a pole and nonpole term, summarized in Refs. [|i"U23f. We only note here that the pole 
term determines the coupling between the NN and ttNN channels. Due to this coupling, 
pion scattering can also be presented as pion absorption on one nucleon and subsequent 
emission from another. Such a two-body mechanism is not included in the standard multiple 
scattering theory. At large momentum transfers this process could be important because of 



the possibility of momentum sharing. However, as argued by Jennings p4| , such a two-body 
mechanism is almost completely canceled by additional contributions which come from the 
time ordering of the pion absorption and emission processes. 

On the other hand, the issue concerning the so-called Pn-problem cannot be regarded 
as settled since the strong cancellation takes place only if there is no interaction between 
the two nucleons in the intermediate state. Our simple estimates show that if we would 
include the repulsive part of the iViV-interaction which involves heavy meson exchange (as 
was done by Mizutani and Koltun in their study of the nd scattering length p5|), the 
cancellation disappears. After a modification of the free two-nucleon propagator by adding 
the mean value of the residual iVTV-interaction < Vnn(p,u, ...) > the corresponding two- 
body contribution becomes proportional to < Vnn >■ Fig. 2 illustrates the effects which 
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could come from this two-body mechanism with < Vnn >= 40 MeV. Clearly, a more 
microscopic calculation is called for, here we want to merely point out that this mechanism 
can improve the situation at backward angles. 



B. Polarization observables. 

Here we begin our discussion with the results for the tensor analyzing powers. In order 
to understand the general structure of these observables and their behavior in different 
kinematical regions we first consider the size of the observables at backward angles 9 W = 180°. 
It is instructive to also neglect the contribution of the deuteron 3 _Di-component. As an 
example, Fig. 3 presents T 20 , r 2 i and r 22 at = 294 MeV. We remind the reader that, 
according to Eq. (19), the last two observables are linear combinations of T 20 , T 21 and T 22 . 

It follows from Eq. (30) that at 9 W = 180° the amplitudes £>, C and V vanish. Therefore, 
at backward angles the observables T 2 i and T 22 vanish as well, and T 20 reduces to the 
following expression 



T 20 (180°) = 



V2(\A 


2 _ 


1 £ 


I 2 ) 


2 


A\ 


2 + l 


£ 1 


2 



(34) 



From Eq. (30) we can see that in the region where the radial integral /% = V2R^ , the 
contribution from the A amplitude becomes very small. In this region we then have 

T 20 (180°) « -V2 , r 22 (180°) « 2r 21 (180°) « — -= . (35) 

v3 

This result is confirmed experimentally for T n > 200 MeV. Note, that at lower energies in the 
region where the deuteron D-w&ve contribution is small (Rqq » R^) tensor observables at 
backward angles become proportional to the strength of the deuteron D-state. For example, 
T 20 (180°) » -4flgV^- 

Another simplified scenario is given by neglecting the pion rescattering contributions 
and the deuteron D-wave. In this case, the C-amplitude vanishes, and ^4^^ = and 
£>( s ) = —T)( s \ Due to these identities, T 2 P = and the other observables can be expressed 
solely via the elementary nN scattering amplitudes 
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2^2 


A s 


\ 2 sm 2 9 w 


12 


A 


2 +8 | 


A s 


2 sin 2 9 n 



T (S) _ ^ I As r sin ^ (5 ) _ 6 ( 5 ) , , 

J 20 i o I /t 12 , o I /I 12 „:„2/i ' J 22 ~\ o J 20 ■ W 



As we have shown in the analysis of the differential cross sections, due to the dominance of 
the p-wave contribution in ttN scattering, Ao(90°) ~ 0. Therefore, in the 9 n = 90° region 
where the contribution from the deuteron S'-state dominates we have 

T 2 ( 5) (90°) « -V2/4, r 2 ( f(90°) w 2r^(90°) w ~. (37) 

Fig. 3 shows that these predictions can already reproduce the experimental measurements 
fairly well. 

The next three figures, Fig. 4, Fig. 5 and Fig. 6, present tensor observables at 7^=134, 
180, 219 and 256 MeV. In general, their behavior follows the simple picture described above. 
At backward angles T 2 o and T 22 are almost completely given by the 3 Di-component of the 
deuteron wave function. In forward direction (up to 9 n = 100°) the main contribution comes 
from the 3 Si-component. Here T 2 o and T 22 contain very little nuclear structure information. 
T21 depends entirely on the deuteron D-state. In the case of the T2i-observable, which 
is mainly associated with T 21 , the deuteron S'-state tends to fill in the minimum around 
Ott = 90°. Note that this minimum appears because of the p-wave dominance of the ttN 
scattering amplitude. Finally, we see that the influence of pion rescattering effects on all 
tensor analyzing powers is small. 

Thus, we can obtain a good description for the tensor observables. This may be related 
to the fact that they are determined mainly by the absolute values of the A — £ amplitudes. 
In contrast, the vector analyzing power iTn depends on the interference between these 
amplitudes. Therefore, we expect this quantity to be more sensitive to details of the model 
ingredients. In Fig. 7 we compare our results with some old |33[ and more recent [32 



measurements. We obtain satisfactory agreement with these data only at T 7r =100 MeV. 
At all other energies our results fail to reproduce the experimental data, especially for 
9^ > 90°. The failure becomes more dramatic with increasing pion energy. In our full 
calculation %T\\ goes through zero around 7^=180 MeV and becomes negative at 9 n > 90°. 
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We have found that pion rescattering is mainly responsible for this effect. Experimental 
measurements show quite the opposite behavior: Above T 7r =250 MeV, iTu changes the sign 
in forward direction producing a negative dip around 6 n = 70° and it assumes large positive 
values for 9„ > 90°. This behavior could only be described in simple calculations without 
pion rescattering (dashed curves). Note that the influence of the deuteron D-state on iT u 
is small. Therefore, in the plane wave approximation the vector analyzing power can be 
expressed via the elementary 71N scattering amplitudes 

■rpiS) _ VElm{A* s A ) sing,, 

11 ~ 12\A \*+8\A s f sin 2 £/ W 

In this paper we have not attempted to present a detailed comparison with the Faddeev 
calculations that have been performed. A review of these results can be found in Ref . |IJ . We 
only mention that all conventional calculations encounter the difficulties in the description 
of the measured vector analyzing power in the A-resonance region. 

On a phenomenological level we have studied the effect of modifying the spin-flip ampli- 
tude As (similar to A at low energies) 

< d I e wK/2£i^: I d> 

A s -> A s - B s => — r - , 39 

<d\ e^-^l 2 \d> 1 ^ 

where | d > is the deuteron ground state. The phenomenological term with a complex 
parameter Bs would be associated with second-order contributions. The sensitivity of the 
differential cross section and tensor analyzing powers to this correction is small since they are 
dominated by the coherent contribution from the non-spin flip amplitude A . However, it has 
a dramatic influence on the vector analyzing power. Using this sensitivity we extracted the 
Bs parameter from the experimental data for iTu. The result of our fit is shown in Fig. 7 by 
the dash-dotted curves. It is remarkable that the energy behavior of the real and imaginary 
parts of Bs follows a typical resonance structure (see Fig. 8). This might be an indication 
for the presence of a residual AN interaction. We point out that a similar situation occurs 
in pion scattering on 3 He for the target asymmetry A y = \[2iT\\ EM. It seems reasonable 



to assume that the origin of this phenomenon is similar for both reactions. It therefore 
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becomes very important to find a consistent description of polarization observables in pion 
scattering on the deuteron and 3 He. 

IV. CONCLUSION 

In this paper, we have studied the interaction of pions with the deuteron in a multiple 
scattering approach carried out in momentum space. Our investigation covered the energy 
region of T n = 60 - 300 MeV, thus covering both the low-energy and the A region. Paris 
wave functions were employed to describe the deuteron. Phase-shift parameterizations were 
used for the elementary irN amplitudes, along with a separable potential for the off-shell 
extrapolation. The full spin and isospin dependence of the 7riV amplitudes were taken into 
account. 

Our framework of multiple scattering theory is clearly less sophisticated than the many 
Faddeev calculations that have studied pion deuteron scattering over the years. The primary 
purpose of our investigation was to demonstrate that the approach developed recently for 
the description of the pion interaction with A > 3 nuclei can be reliably applied in the case 
of the deuteron. The results obtained are transparent and can in the following be employed 
for the study of pion photo- and electroproduction off deuterium. 

We found that our approach gives a good description of the differential cross sections 
and the tensor analyzing powers. At low pion energies, second-order corrections must be 
included, while at higher energies and backward angles discrepancies appear that may be 
related to two-body processes involving heavy meson exchange. The tensor analyzing powers 
strongly depend on the deuteron D-state and can be reproduced satisfactorily. Clearly, the 
outstanding problem is the explanation of the well-measured vector analyzing power which 
has defied a number of attempts using more advanced approaches than ours. While there 
may be a sensitivity regarding a residual AiV-interaction care has to be taken to include it 
consistently. 

In the following paper (Part II), we will show that the pion-deuteron interaction devel- 
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oped here is adequate to describe the final state interaction in coherent pion photoproduction 
from deuterium. 
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In this Appendix, we give the complete expressions for the Robson amplitudes A,...,S 
where the full contributions from the deuteron D-state have been taken into account. 

First, we express the nuclear matrix elements M SLJ (Q) defined by Eqs. (25-26) via the 
radial integrals of Eq. (27) 



APPENDIX: 




(Al) 




(A2) 




(A3) 




(A4) 



Then the amplitudes A, ...,£ in the PWIA approach can be written as 



A=\ —A W Mooo(Q)-— 7 =(l-3cos^)M 022 (g) 




(A5) 
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B = - 



1 4tt 

2YT smfl " 



A 5 W s M S (Q) - - A W s M 022 (Q) 



C = -- J — (1 + cos^) A M 022 (Q) , 



^ 1 / 47r • n 

v = 2\T Bm6ir 



A s W s M S (Q) + - A W s M 022 (Q) 



MoooW) + ^=(1 - 3cos^)M 022 (Q) 



In Eqs. (A6) and (A8) we introduced the matrix element 



M S (Q) = M W1 (Q) - -j=M 121 (Q) . 
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FIGURES 

FIG. 1. Differential cross sections for tt + elastic scattering on the deuteron at pion kinetic 
energies T n =65, 181 and 254 MeV calculated with (solid curves) and without (dashed curve) spin 
flip transition. The dotted curves at 7^=181 and 254 MeV are the results obtained without the 
deuteron D-state. The dash-dotted curve at T n =65 MeV is the result with second order corrections 
(32). Experimental data are from Refs. §6|,|7|(») and Ref. @(o). 



FIG. 2. Differential cross sections for ir + elastic scattering on the deuteron at pion kinetic 
energies 7^=142-324 MeV. Solid and dashed curves are full and PWIA calculations, respectively. 
The dotted curves are our best fit with a phenomenological interaction of two intermediate nucleons 
in the Jennings mechanism. Experimental data are from Ref. ^7)(»), Ref. f28|(o) and Ref. j29|(A). 

FIG. 3. Tensor analyzing powers T20, T21 = T21 + |(722 + T20/V6) and T22 = T22 + 
at T 7r =294 MeV (solid curves). Individual contributions from T21 (in T21) and T22 (in T22) are 
shown by dashed curves. The dotted curves are the results obtained for T20, 721 and T22 without 
the deuteron 79-state. Experimental data are from Ref. |3(|(»), Ref. fl31|| (black triangles) and Ref. 
Wo). 



FIG. 4. T 20 observable at 7^=134-256 MeV. Solid and dashed curves are full and PWIA 
calculations, respectively. The dotted curves are the results obtained without the deuteron D-state. 
Experimental data are from Ref. (33l](») and Ref. p0|(o). 



FIG. 5. The same as in Fig. 4 for the observable T21 = T21 + \ (Ti2 + 7 , 2o/v / 6)- Experimental 
data are from Ref. @|(») and Ref. |30[|(o). 



FIG. 6. The same as in Fig. 4 for the observable T22 = T22 + T20/VQ. Experimental data are 
from Ref. fl)(«) and Ref. ||(o). 
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FIG. 7. Vector analyzing power iT n at a) 7^=100-164 MeV b) 7^=180-294 MeV. Solid and 
dashed curves are full and PWIA calculations, respectively. Dash-dotted curve are the results of 
our fit with phenomenological term (39). Experimental data are from Ref. [^] (•) and Ref. |32|| (o). 

FIG. 8. Energy dependence of the modification for the spin-flip interaction (see Eq. (39)) 
extracted from the vector analyzing power iTn. 
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